Anderson localization from the replica formalism 
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We study Anderson localization in quasi-one-dimensional disordered wires within the framework 
of the replica <r-model. Applying a semiclassical approach (geodesic action plus Gaussian fluctua- 
tions) recently introduced within the context of supersymmetry by Lamacraft, Simons and Zirnbauer 
Pl, we compute the exact density of transmission matrix eigenvalues of superconducting wires (of 
symmetry class CI.) For the unitary class of metallic systems (class A) we are able to obtain the 
density function, save for its large transmission tail. 

PACS numbers: 72.15.Rn 



At present, there exist two theoretical approaches ca- 
pable of describing strongly localized phases of disordered 
wires: supersymmetry_(SUSY) [2j and the DMPK trans- 
fer matrix approach This represents a serious limi- 
tation in as much as both formalisms are ill-suited for 
generalization to the presence of Coulomb interactions 
(see, however, Ref.Q.) Reciprocally, it has, so far, not 
been possible to describe strong localization phenomena 
by those theories that may be applied to the analysis 
of interaction effects — replica field theory [j| and the 
Keldysh approach 0. 

It is the purpose of this letter to introduce a replica 
field theory approach, capable of describing strongly lo- 
calized phases. Conceptually, our work is based on a 
recent paper yij by Lamacraft, Simons and Zirnbauer 
(LSZ) where saddle-point techniques have been applied 
to analyze the SUSY generating functionals of quasi one- 
dimensional disordered conductors. Specifically it was 
shown that four out of ten symmetry classes of disor- 
dered metals are semiclassically exact pf in that the sta- 
tionary phase results coincide with those obtained by 
DMPK methods 8| . We here show that the phenomenon 
of semiclassical exactness pertains to the replica formal- 
ism and, in particular, 'survives' the analytical continu- 
ation inherent to that approach. Applying the technique 
to the non-semiclassically exact unitary symmetry class, 
we find that it still produces qualitatively correct results. 

To introduce the replica-generalization of the method 
we consider a disordered superconducting wire in the 
presence of spin-rotation and time reversal invariancc 
(symmetry class CI in the classification of Ref. 9J.) The 
(thermal) transport properties of this system may be con- 
veniently characterized in terms of the average density 
of transmission matrix eigenvalues, p((f>). Within the 
fermion-replica formalism the latter may be expressed 
through the generating function 

R 

Z(0) = []det(l-sin 2 (0 a /2)ttt), 

a=l 

where tt^ is the transmission matrix with eigenvalues 



Tj = cosh~ 2 (0j/2) and 9 = diag(#i, . . . , Br). Defining 
the function F(9) = lim^^o^| e _ | ^(^), the trans- 
mission matrix eigenvalue density is obtained as |lfj| : 
p(<j>) = ±(F(i<t> + ir)-F(i<fi-n)). 

The field theoretical representation of the generating 
function for class CI is given by 

g{T) t 
Z(9) = jvg e- s ^\ S[g] = ~ Jdt tv{dgdg- 1 ) , (1) 

9(0) 

where g is a field of matrices g(t) € Sp(2i?), the func- 
tional integration extends over the Haar measure on the 
symplectic group, and T = L/t; is the length of the wire, 
L, in units of the localization length £. At the left and 
right end point of the wire the field is subject to bound- 
ary conditions flfl lll] | which in the case of class CI read 
as <?(0) = 1 and g(T) — exp(i9(g) (73). Here the Pauli ma- 
trix 03 acts in the space defining the symplectic condition 

Our strategy will be to subject the functional to a 
straightforward stationary phase analysis |12| . Varying 
the action S[g] w.r.t. g, one obtains the Euler-Lagrange 
equation: S g \ g= -S[g] = => d(g~ 1 dg) = 0, which inte- 
grates to the condition g~ 1 dg = const. The solutions to 
this latter equation are given by g — exp(iWt/T), with 
constant Lie-algebra elements W <E sp(2R). Evaluating 
g at the system boundary t = T, we obtain the condition 
exp(iW) = exp(i9(g>a 3 ). This is solved by W = 9^®az, 
where 9^ = 9 + 27rn and h = diag(rti, ... ,tir) is a vec- 
tor of integer 'winding numbers'. The saddle point action 
is given by %W] = ^ ^f =1 (6»i n) ) 2 , indicating that at 
length scales, T > 1, mean field configurations travers- 
ing multiply around the group manifold become energeti- 
cally affordable. Physically, these configurations describe 
the massive (and perturbatively inaccessible) buildup of 
interfering superconductor diffusion modes. Their pro- 
liferation at large length scales forms the basis of the 
localization phenomenon. 

To obtain the contributions of individual saddle points, 
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g^ n \ to the generating function, we need to integrate 
over quadratic fluctuations. We thus generalize to field 
configurations g(t) = exp(iW(t)) jW, where the fields 
W(t) <E sp(2i?) obey vanishing (Dirichlct) boundary con- 
ditions W(0) = W(T) = 0. Parameterizing these fields 

as W = Su = oW/i ® °>> where cto = I2 and W 7 ^ are 
R x R hcrmitian matrices subject to the Lie algebra 
constraints W = -W T and If, = Wf, i = 1,2,3, 
the quadratic expansion of the action reads as: S[g] — 
% (n) ] + Si[W , W 3 ] + S II [W 1 ,W 2 ] + 0(W 3 ), where 



Si[W ,W 3 ] = ~ [dttr(dW dW + dW 3 dW 3 
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{W dW 3 + W 3 dW ) §W Y 



S/j[Wi,W 2 ] = ^ JdttY(dW 1 dW 2 + i e -^-W l dW 9 ( ' : 



The integration over the matrices W M leads to fluctuation 
determinants, which may be calculated by the auxiliary 
identity det(-<9 t 2 + 2zT~ 1 d t ) = smh(z)/z, where z G C, 
and the differential operator acts in the space of functions 
obeying Dirichlct boundary conditions. As a result we 
obtain the stationary phase generating function 
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where the first/second fluctuation factor stems from the 
integration over the field-doublets (Wo, W 3 )/(Wi, W2). 
(In passing, we note that as an alternative to 
the brute force integration outlined above the re- 
sult (J2J) can be obtained by group theoretical rea- 
soning: according to general principles [lif. the fluc- 
tuation integral around extremal (geodesic) configura- 
tions on a general semi-simple Lie group is given 
by : ria>o«(ff ( " ) ) sin " 1 («(5 ( " ) ))exp(-5[5(™)]), where 
the product extends over the system of positive roots 
of the group, a(g^). Equation J5J above is but the 
Sp(2i?)-variant of this formula.) 

In the limit of coinciding boundary phases, 9 a — > 9, 
the denominators sin[(#i n ' ) — 9^ )/2\ — > 0, i.e. the con- 
tribution of configurations h containing non-vanishing 
winding number differences n a — n a ' 7^ diverges. (At 
the same time, we do know that the integration over the 
full group manifold must generate a finite result. In- 
deed, it turns out that if we first sum over all winding 
number configurations n and only then take the limit of 
coinciding phases, all divergent factors disappear.) This 
divergence reflects the presence of a zero mode in the 
system: for uniform boundary phases, 9 cx transfor- 
mations — > exp(iV°)g^ exp(— iV°) with constant 
block-diagonal V° = V® £g> Co + U 3 ° <£> 03 conform with the 
boundary conditions but do not alter the action. 

As we shall see below, the presence of zero modes im- 
plies that only winding number configurations of the spe- 
cial form (n, 0, . . . , 0) survive the replica limit, R — ► 0. 
However, before elaborating on this point, let us eval- 
uate the contribution Z n of the distinguished config- 
urations to the generating function. Throughout we 
will denote the boundary angles by 9 a = 9 + rj a , un- 



derstanding that the limit 77^ — > is to be taken at 
some stage. (Within this representation, the 'free en- 
ergy' F(9)d 9l \ 9 ^ e Z(9) = d m \ v ^ Q Z(9 + r)).) The 'dan- 
gerous' product [] a<a ,(...) in Eq. then reduces to 
~ (7rn/sin(?7i/2))- R_:L » (2Trn/r]i) R ~ 1 ; all other contri- 
butions to Z n are finite. The appearance of a pole of 
(R — l)st order hints at the presence of R — 1 complex 
zero modes (generated by the R — 1 components of the 
matrix V° that do not commute with g( n ' .) At this stage, 
we take the limit R — > 0. As a result, the divergent fac- 
tor gets replaced by a 'pole of degree (—1)', i.e. the zero: 
r)i/(2irn). (It is worth noting that in SUSY a contribu- 
tion similar to the singularity of degree (—1) is obtained 
by integration over the non-compact bosonic degrees of 
freedom; the complementary single replica channel a = 1 
corresponds to the fcrmionic sector.) Therefore the sub- 



sequent differentiation (F[ 
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+0 -Z) must act on 



this linear factor 771 , all other occurrences of rj a in Z n 
may be ignored. 

Evaluating the partition function in this manner, we 
obtain Z„ #0 = ^r^^exp(-7rn(7rn + 9)/T). We fi- 
nally differentiate w.r.t. r\\ and arrive at the result 
p(<p) = Po(4>) +J2n=io Pnify) where the 'Drude plus weak 
localization term' p = (2T) _1 - ((j> 2 + Tr 2 y 1 /2, while the 
non-perturbative contributions are given by: 



p— Tp-n(n-l-l) 

Pn{<l>) = - —^-o Rc 

27T Z 71 



+ i7r(2n + l) j2 n£ 



ip + iir(n + 1) 



(3) 



This expression identically coincides with the SUSY re- 
sult 0, and with the exact DMPK result @|. To illustrate 
the 'crystallization' of the transmission matrix eigenval- 
ues at the discrete values 4>j ~ 2j'T, the function p(cj)) 
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FIG. 1: Top: Density of transmission eigenvalues of the super- 
conductor class CI for T=.02 (dashed-dotted); 1 (dashed); 50 
(full). Bottom: the same for the unitary class A. The negative 
density at small <f> represents an artefact of the saddle-point 
approximation. 

is plotted in Fig. for a few values of T. Following 
LSZ, the heat conductance of the wire may be obtained 
by integrating the result above against the weight func- 
tion 1/ cosh 2 (0/2). Summing the result of this integra- 
tion over winding numbers, one obtains the asymptotic 

result 1] g T ~ 1 Ae- T /V^T . 

Our so far analysis focused on the specific set of wind- 
ing number configurations, (n, 0...,0). To understand 
why contributions of different structure vanish — a fact 
that greatly simplifies the formalism — consider the set 
(0, . . . , n, . . . , 0). By symmetry, winding number config- 
urations of this type will lead to an expression similar 
to Z n above, only that the leading pre-factor gets re- 
placed: 77i/(27rn) — > r) a /(2irri), where a £ {2, . . . , R} 
marks the position of the non-vanishing winding num- 
ber. Since, however, we still differentiate w.r.t. rji, 
this contribution vanishes in the limit rj a — > 0. The 
argument above may be generalized to generic contri- 
butions, (n%, ri2, ■ ■ ■ , tir) (n, 0, . . . , 0). (By symmetry, 
one may order the winding numbers in an ascending or- 
der (0, . . . , 0, 1, . . . , 1, 2 . . .). Assuming that there are N n 
winding numbers n (where N n — R) and choosing 
the boundary angle in the sector n to be 9 + nrj, one 
verifies that for any fixed configuration, the R — > re- 
sult contains uncompensated powers of rj and, therefore, 
vanishes.) 

Before proceeding, it is worthwhile to compare the 
mean field analysis above to the more established field 
theory transfer matrix method 0. To this end, let us 



interpret Z{9) = (g(T)\exp(-TH)\t) as the path inte- 
gral describing the (imaginary time) quantum mechani- 

cal transition amplitude |1) — ► \g(T)) of a particle on the 
group space Sp(2i?). The Hamiltonian corresponding to 
the (purely 'kinetic') action of the path integral is given 
by H = — 2 A where A is the Laplace operator of the 
group space Sp(2i?). 

Our analysis above has been tantamount to a semiclas- 
sical or WKB analysis of the transition amplitude. Alter- 
natively, and more rigorously, one may employ the spec- 
tral decomposition, Z{9) = Ylx^xid^^W exp(—Te\), 
where tp\ are the eigenfunctions of the Laplace opera- 
tor, e\ its discrete energy eigenvalues and g = g(T). 
For general Lie groups (and supergroups) formal expres- 
sions for these spectral decompositions are known [l4| . 
Noting that for large systems I > ^, only eigenstates 
with minimal energy e\ effectively contribute to the sum, 
this knowledge has been used to compute the localization 
properties of disordered quantum wires within the SUSY 
formalism 0, EI • The problems with transferring this 
approach to the replica formalism lie with the analytical 
continuation from integer group dimension R to R — ► 0. 
In taking this limit, it is essential to keep track of high- 
lying contributions to the spectral sum. These terms 
grow rapidly more complex which is why attempts to ob- 
tain a replica variant of the 'quantum approach' above 
have failed so far. 

Having discussed the method for a symmetry class 
that enjoys the semiclassical exactness, we next out- 
line what happens in cases where this feature is ab- 
sent. By way of example, consider a metallic disor- 
dered quantum wire in the absence of time-reversal in- 
variance — the unitary symmetry class, A. In this 
case, the fermionic replica generating function is given 

by Z{9) = VQ exp (-\ ^dttv{dQf\ where the 

matrix Q(t) E U(2i?)/U(i?) x U(R) H 0, and the 
boundary configurations are given by Q(0) — 03 <S) 1 and 
Q(T) = e- iCT2 ® e "/ 2 cr 3 e ia2 ® § / 2 . Here, the two-component 
structure distinguishes between advanced and retarded 
indices. As before, the stationary phase configurations: 
Q(t) = e -i°2®^ n) t/(2T) a3 e i CT2 ®^>t/(2T) of the functional 

integral do not mix different replica channels. Geometri- 
cally, they can be interpreted as trajectories (in general, 
with non-zero winding number, n) on the meridian of the 
sphere U(2)/U(l) x U(l) (the single replica manifold.) 
Fluctuations may be conveniently parameterized by gen- 
eralization cr 3 — > e m ^a- ln where W = W\ <8> 01 + W2® a% 
and W\ t 2 are hermitian R x R matrices. 

The subsequent calculations largely parallel those for 
class CI above. Expanding to second order in W\ t 2 and 
performing the Gaussian integration, we again observe 
that only winding number configurations (n, 0, . . . , 0) 
survive the analytical continuation procedure, R — > 0. 
Differentiating w.r.t. 6\ and then putting 6 a — > 9, we 
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obtain the result: p{cj)){2T)- 1 ~ £ n#0 (-l) n p„(</>), where 



2n 2 n 



Re 



y/(4>+iir)(cf)+iir(2n+l)) , w 
<f> + m(n + 1) 



(The same result is obtained by saddle-point analysis of 
the SUSY generating functional.) In Fig. ^p, the func- 
tion p{4>) is plotted for several values of the parameter T. 
For small T the density is almost constant, reflecting the 
Dorokhov distribution of eigenvalues 0, ^| . For large 
values of T the spectrum crystallizes at <j>j « (1 + 2j)T. 
The lowest eigenvalue </>o does, indeed, correctly deter- 
mine the localization length of the system. Except for 
the evident failure of the method at small values 4> <C (f>o 
[lif . the large scale profile of the DoS is in good agree- 
ment with results obtained by the transfer matrix meth- 
ods munis. 

Summarizing we have shown how the localization phe- 
nomenon in quasi one-dimensional systems may be de- 
scribed by a semiclassical approach to fermionic-replica 
field theories. We were able to reproduce the exact trans- 
mission matrix eigenvalue density for symmetry class CI, 
while for the unitary class we obtain qualitatively correct 
results (except for the tails of the eigenvalue spectrum.) 
The comparative simplicity of the approach makes us be- 
lieve that it may be successfully applied to problems that 
can not be treated by other means. Evidently, the next 
direction of research will be the study of the impact of 
Coulomb interactions on the localization phenomenon. 
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